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Abstract 

Transgenerational effects are broader than only parental relationships. Despite 
mounting evidence that multigenerational effects alter phenotypic and life-his- 
tory traits, our understanding of how they combine to determine fitness is not 
well developed because of the added complexity necessary to study them. Here, 
we derive a quantitative genetic model of adaptation to an extraordinary new 
environment by an additive genetic component, phenotypic plasticity, maternal 
and grandmaternal effects. We show how, at equilibrium, negative maternal 
and negative grandmaternal effects maximize expected population mean fitness. 
We define negative transgenerational effects as those that have a negative effect 
on trait expression in the subsequent generation, that is, they slow, or poten- 
tially reverse, the expected evolutionary dynamic. When maternal effects are 
positive, negative grandmaternal effects are preferred. As expected under Men- 
delian inheritance, the grandmaternal effects have a lower impact on fitness 
than the maternal effects, but this dual inheritance model predicts a more com- 
plex relationship between maternal and grandmaternal effects to constrain phe- 
notypic variance and so maximize expected population mean fitness in the 
offspring. 
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Introduction 

Maternal effects occur when maternal phenotype(s) influ- 
ence offspring phenotype(s) by means other than direct 
genetic transmission (Mousseau and Fox 1998). This 
alternative mode of inheritance can act as a divergent or 
stabilizing force (Rasanen and Kruuk 2007; McGlothlin 
and Galloway 2014) by, respectively, accelerating adapta- 
tion to novel environments (Lande and Price 1989) or by 
reducing phenotypic variance to maximize population 
mean fitness in relatively stable environments (Hoyle and 
Ezard 2012). When the environment experienced by the 
parent covaries with the environment encountered by the 
offspring, maternal effects are particularly likely to evolve 
(UUer 2008; Fischer et al. 2011; Kuijper and Johnstone 
2013). Transgenerational effects are not restricted to the 



maternal generation, however. Phenotypic "memory", 
often provoked by environmental factors, can persist for 
many successive generations (Molinier et al. 2006) and 
may differ depending on the sex of the influencing ances- 
tor (Lock 2012). In particular, levels of grandmaternal 
nutrition can, over and above maternal effects, alter life- 
history traits in the coUembolan Folsomia Candida (Hafer 
et al. 2011), obesity in humans (Cropley et al. 2006), and 
mass in white-tailed deer Odocoileus virginianus (Mech 
et al. 1991). Analogously, a grandfather effect was 
reported for offspring weight in the western mosquitofish 
Gambusia affinis (Reznick 1981). The food environment 
in the great-grandmaternal generation in soil mites San- 
cassania berlesei can leave a phenotypic signature on life- 
history traits including egg length (Plaistow et al. 2006). 
WhUe the potential for multigenerational transgeneration- 
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Maternal & Grandmaternal Indirect Genetic Effects 

al effects on life-history and phenotypic traits is therefore 
dear, it is not known how these effects from different 
generations should combine to maximize population 
mean fitness in the offspring generation. Here, we extend 
recent modeling work (Lande 2009; Hoyle and Ezard 
2012) to find out. 

We assess fitness within a focal generation, whose pheno- 
types are constructed from multiple components. The 
transgenerational maternal and grandmaternal components 
assume that the phenotypes in previous generations deter- 
mine part of the current phenotype under selection and so 
the fitness belongs (Wolf and Wade 2001) to the current 
generation. Wolf and Wade (2009) defined a "true" mater- 
nal effect as one with a causal link between maternal geno- 
type or phenotype and the offspring phenotype. Within 
this definition, various subcategories have emerged: Trans- 
generational effects that enhance fitness have been termed 
"adaptive maternal effects" (Marshall and UUer 2007), 
while Rasiinen and Kruuk (2007) defined positive maternal 
effects as ones that (could) accelerate microevolution, pre- 
sumably toward some shifting optimum. The challenge 
with more specific terms is that their combinations are not 
fixed: Hoyle and Ezard (2012) showed how a negative 
maternal effect is "adaptive" in the sense of Marshall and 
UUer (2007) in a relatively stable environment, but a posi- 
tive maternal effect is "adaptive" during adaptation follow- 
ing a rapid shift in environment. We therefore do not 
assume that a particular maternal effect is adaptive per se 
because this definition requires an understanding of the 
context in which the organism lives (Plaistow et al. 2006; 
Plaistow and Benton 2009). Our goal is a better under- 
standing of how positive and negative maternal and grand- 
maternal effects result in changes to expected mean fitness 
in the offspring. 

The basic framework is a quantitative genetic model of 
adaptation via plasticity, genetic assimilation (Lande 
2009), maternal (Kirkpatrick and Lande 1989; Hoyle and 
Ezard 2012), and grandmaternal effects. Lande (2009) 
used this approach to study how phenotypic plasticity 
and additive genetic variance interact during adaptation 
to an extraordinary new environment. Hoyle and Ezard 
(2012) blended this model with that of Kirkpatrick and 
Lande (1989) to include a maternal effect coefficient m 
that characterized the direct path from maternal pheno- 
type to offspring phenotype, separately from the genetic 
contribution. Note that this maternal effect is therefore 
an indirect genetic effect because the maternal phenotype 
is, in part, genetically determined and the inheritance 
pathway is not via direct genetic transmission (Moore 
et al. 1997; McGlothlin and Brodie 2009). We extend the 
model in Hoyle and Ezard (2012) to include a linear con- 
tribution of the grandmaternal phenotype, captured via 
an analogous (indirect genetic) grandmaternal effect coef- 
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ficient g. Note that we model an explicit grandmaternal 
effect, as demonstrated experimentally (Hafer et al. 2011; 
Lock 2012), and not the cumulative consequences of suc- 
cessive maternal effects. We investigate how m and g 
combine to determine fitness in two environmental sce- 
narios: at equilibrium and in the immediate aftermath of 
a sudden environmental shift, because a traditional per- 
spective of maternal effects is that they facilitate adapta- 
tion to a changing environment (Rasanen and Kruuk 
2007; UUer 2008). 

Materials and Methods 

We take a quantitative genetic approach. The model is 
univariate, asking how a trait affects the same trait in 
future generations (e.g.. Falconer 1965). Generations are 
discrete and nonoverlapping. We assume phenotypes 
determine fitness and are constructed from multiple com- 
ponents. The phenotype of an individual at time (=gen- 
eration) t is: 

z, = at + fot£t-T + rnz*_^ + gz*_^ + e,, (1) 

where a, is the additive genetic component (breeding value) 
of the phenotype in the reference environment (e = 0) and 
bf is the linear slope (reaction norm) of the plastic pheno- 
typic response to the environment (Lande 2009). We note 
that the mean breeding value is, by definition, constrained 
to have a mean of zero. There is a lag t, measured in frac- 
tions of a generation, between juvenile development and 
the time when selection occurs. This lag is a key parameter 
in disentangling within-generation phenotypic and trans- 
generational plasticity (UUer 2008; Hoyle and Ezard 2012). 
The maternal coefficient m and the grandmaternal coeffi- 
cient g characterize the influence of the maternal and 
grandmaternal phenotypes after selection (z*_j and z*_2, 
respectively), on the offspring phenotype. is the residual 
component of phenotypic variation, which we assume to 
have mean zero without loss of generality. Fig. 1 is a sche- 
matic of model structure. 

From equation 1, the variance of Zf in a constant envi- 
ronment £ satisfies: 

ffz, ~ G„a + GhbE^ + 2mGa,zi^ + '^SGa,zi^ + 2mGi,,z;_^e 
+ 2gGh,zi,e + m^(y\> ^ + IrngG^i^zi^ + g^^l-^ ^ + '^e' 

where Gaa and Gyi, are the variances of a, and b,, respec- 
tively, and are assumed to be constant. Ga,z* is the 
covariance of a, with z*_^; aU other covariances are 
defined in simUar fashion. Following Lande (2009), the 
phenotypic variance is minimized in the reference envi- 
ronment £ = 0, which forces the covariance between addi- 
tive genetic effects and phenotypic plasticity Gab to be 
zero. 
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Figure 1. Schematic of model structure showing how Z(, the 
offspring phenotype in generation t, is constructed from additive 
genetic at and phenotypically plastic bt components. The 
transgenerational effects are modeled using the maternal coefficient 
m and an analogous grandmaternal coefficient g, which link the 
maternal and grandmaternal phenotypes after selection {Zt_, and 
Zf_2, respectively) to Zf. et-z the environment before selection. 

Also, from equation 1, we see that 

Ga,z, =Gaa + »jGa,z* , +^Go,z«_^, 

The mean phenotype (an overbar denotes the expecta- 
tion value) is given by 

z, = a, + btEt-z + mz*_i + gz*_^, 
Finally we use 

z* = a* + h*e.+ mz*_^^ + gz*^^ + 

where * denotes values after selection in generation t, to 
see that 

Gz,-_,z,-_, = G„;__z,-_, + Gfc;_,z,-_^,£+ +^Gz;_^z,-_3 

Now, we look for an equilibrium under weak selection, 
such that: 



4Ga,z;_ 
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2Ga z* , — 2Ga; 
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(7. = a^. 



We assume that at equilibrium, under weak selection, 
the distribution of e* is the same as that of and hence 

G.. ,z- _ = 0. 
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The optimum phenotype is assumed to be a linear 
function of the environment: d, = A + Be, where A and B 
are constants and e, is the environment at time t. The 
fitness is defined as a Gaussian function given by 



W{£,,Zt) = Wn,ax CXp |- ^ ' | 



where W^^^ and oj are constants. Here, oj^ is the 
strength of stabilizing selection (width of the fitness 
function). 

Assuming that the phenotype Zt follows a Gaussian 
distribution, the population mean fitness is given by 



W(et,z,) = 



yffl2exp|-^(zf-0f)^| 



where v = 



Following Lande (1979), and using 
z*_j = Ut + bt£t-i-T + ntz*_2 + gz*_^ (which assumes no 



fertility selection), we therefore have 

Aa = -yGaa{at -A + b,£,_^ - Be, + mz*_^ + gz*_2) 
(1 + m), 

Ab = -yGbbia, - A + btSt-z - Bs, + mz*_^ + gz*_2) 



(2) 



(3) 



(et_i + m£t_i_0, 
z* = at + Aa + {b, + Ab)e,_, + mz*__j + gz*_2. (4) 

Adaptation to new environment 

We consider —0.6 < m < 0.6 and —0.3 < ^ < 0.3 at incre- 
ments of 0.05 to investigate how m and g interact to 
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determine fitness after a sudden environmental shift, 
modeled as a noisy step change e, = U,d + (Lande 
2009). In this model of environmental change, [7, is the 
unit step function that jumps from 0 to 1 at f = 0, (5 is 
the size of the sudden change in average environment, 
and is a Gaussian stationary autocorrelated random 
process with mean zero, variance cr?, and autocorrelation 
Pt: over the interval t. 

Taking the expectation of equations (2-4) over the 
Gaussian distribution of the stochastic component in the 
environment we have 

E{Aa) = -yfi„a{a, -A+ 8Ut-,bt - ^UtB + mE{z;_^) 

+ gE{zl,))il + m) (5) 

E{Ab) = -yfihb{at -A + - SUtB + mE{z*_^) 

+ gE{zl^)){8U,-, + m5Ut-i-,) 
- yfiuialiMl + m^) - p,B} (6) 

E{z*) = a, + 5Ut-,bt + m£(z*_i) + g£(z*_2) 

- y,GbbSU,-,cjl{b,{l + n?) - p,B} 

- y,Gbb{SU,-, + m5Ut-i-,){b, - p,B)(j\ 

- )'.{Ga„(l + m) + Gbb{5Ut-,{5U,., + m5U,.i-,) 
+ (7^)} X {a, - A + dU,-,b, - dU,B + mE{z*_^) 

+ gE{zl2)}- (7) 

To get the above equations, we assumed that the sto- 
chastic component of the environment is uncorrelated 
across more than a single generation. We have also 
approximated y by = ^ji^\^^2y and, under weak selec- 
tion, approximated the expected phenotypic variance 
E{al) as: 

^Ga, + GbbS'Ul^ + ffl 

-I- Gfcicr. + (G„a + Gbbo Ut-xUt-T-i) 

- 4 — 2m — ^ 

+ - r^^ (Gaa + GbbS^Ut-.Ut-r-l) 

4 — 2m — g 

+ JZ . ^„ r(Gaa + Gbb^^Ut-T-iUt-T-i). 

{\ - g)[A: - 2m - g) 

(8) 

Note that this is slightly different to Hoyle and Ezard 
(2012), in that we now incorporate distinct environments 
after juvenile development (but before selection) in the 
present, maternal and grandmaternal generations rather 
than just the environment in the present generation. We 
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expect these to be good approximations for a\, a\ < <a/. 
The expected value of fitness is 

E{W)= W^axV7^exp{-|£(zt-e,)'}, 

X exp I - y (flf - A + 5U,.,bt - 8UtB 
+ mE{-zl^) + gE{z;_,)f} 

2 

X exp { - + + ^2) +5^ - 2Mp.)}. 

(9) 

Commented MATLAB simulation routines are available 
as online supplementary material. 

Results 

Multigenerational transgenerational effects change the 
dynamics during both transient and equilibrium phases 
(Fig. 2). The potential for maternal and grandmaternal 
effects to accelerate adaptation to the extraordinary new 
environment is clear: Adaptation is fastest when m > 0 
and g > 0 and slowest when m < 0 and g < 0, with m 
being more influential than g (Figs. 2, 3). As an example, 
it takes 5270 generations for expected mean fitness E{W) 
to be within 0.001 of the equilibrium in the novel envi- 
ronment when m = 0 and g = 0. If m = 0.05 and g = 0, 




0.02-1 

1 1 1 1 

0 100 1000 10000 

Time (= generations) 

Figure 2. Expected mean fitness (as a proportion of M/max = 1) 
tinrough time, plotted on the scale of the natural logarithm. All four 
combinations of m ± 0.35 and g ± 0.15 are depicted to illustrate 
how transgenerational effects combine during adaptation to an 
extraordinary new environment. Black if m > 0 and gray if m < 0; 
solid lines if g > 0 and dashed lines if g < 0. Positive 
transgenerational effects lower expected mean fitness (Hoyle and 
Ezard 2012). Parameter values follow Lande (2009) and Hoyle and 
Ezard (2012): A = 0, 6 = 2, = 0.25, (7^ = 2, G^, = O.S, 
Gtb = 0.045, 7 = 0.02, co^ = 50, and <5 = 10. 
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-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

Maternal effect coefficient 

Figure 3. (A) Expected mean fitness (on the scale of the natural 
logarithm) after 100,000 generations and (B) number of generations 
to recover equilibrium fitness after the step change in the 
environment (within 0.0001; see also Fig. 2). The dotted lines in panel 
(A) indicate maximum expected mean fitness holding either maternal 
or grandmaternal effects constant. When g = 0, the model reduces to 
the results of Hoyle and Ezard (2012). Parameter values as in Fig. 2. 

then this point is reached after 4928 generations. If m = 0 
and g = 0.05, then this point is reached after 5088 genera- 
tions. 100,000 generations after the step change, expected 
mean fitness E{W) is maximized for m = —0.4 and 
g = —0.15 (Fig. 3A). Although m > 0 and ^ > 0 accelerate 
adaptation during the transient phase (Figs. 2, 3B), it 
realizes lower EiW) once the new equilibrium is reached. 
At the new equilibrium (after the step change), m < 0 
and ^ < 0 maximize E(W). When either m or ^ is posi- 
tive, however, the optimal ^ or m (respectively) to maxi- 
mize E{W) (dotted lines in Fig. 3A) is of the opposite 
sign. 

The shape of the plateau on the fitness landscape (Fig. 
3A) suggests two key points to explain these results: (1) 
E{W) is more sensitive to the phenotype in the maternal 
generation than the grandmaternal one, particularly if 
m > 0 (compare the distribution of contours in Fig. 3A) 
and (2) the fitness is largely driven by the combination of 
multigenerational components that minimize the expected 
phenotypic variance E{al) (compare Figs. 3A and 4). The 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 
IVIaternal effect coefficient 



Figure 4. Expected phenotypic variance in the reference (A, s = 0) 
and novel environment (B, e = 10). The dotted lines indicate minimum 
expected phenotypic variance holding either maternal or 
grandmaternal effects constant. Parameter values as in Fig. 2. 

variance in the reference environment is lower than in 
the novel environment, but the combinations of m and g 
that influence it are qualitatively similar (Fig. 4). 

Discussion 

This is the first study to investigate how multigenerational 
indirect genetic effects combine to determine expected 
mean fitness. We do this using a straightforward concep- 
tual extension of recent theoretical quantitative genetic 
models (Lande 2009; Hoyle and Ezard 2012). Hoyle and 
Ezard (2012) showed that negative maternal effects maxi- 
mize population mean fitness in relatively stable environ- 
ments because m < 0 minimizes phenotypic variance, 
which keeps more of the population closer to the target 
phenotype. On an individual level, if m is negative, the 
effect of the maternal phenotype is discounted against the 
inherited genes and so the phenotype is brought closer to 
the optimum. The straightforward extension of this logic 
to multiple generations here realizes nontrivial mathemat- 
ical expressions (equations 7 & 8), suggesting that the 
maternal m and grandmaternal g exert complex effects on 
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the expected phenotypic variance and therefore on pheno- 
typic evolution. In the immediate aftermath of a sudden 
environmental shift, adaptation to the new optimum is 
fastest when both m and g are positive (Figs. 2, 3B). Once 
the new equilibrium is reached, both m and g should be 
negative to maximize expected mean fitness _E(W) (Fig. 
3A). 

Experimental work is increasingly documenting how 
grandmaternal effects alter life-history (Hafer et al. 2011; 
Lock 2012) or phenotypic (Cropley et al. 2006) traits over 
and above maternal effects, but, to our knowledge, experi- 
ments assessing the fitness implications after controlling 
for multigenerational factors have not yet been per- 
formed. As expected under Mendelian inheritance (Galton 
1897), E{'W) appears more sensitive to changes in m than 
g during both the transient and new equilibrium phases. 
This sensitivity is more complex than a half contribution 
from mothers and a quarter contribution from grand- 
mothers, however (equation 8). In particular, at equilib- 
rium, the contributions of the two generations are not 
uniform across all explored combinations of m and g: 
Note the vertical and horizontal spread of contours in 
Fig. 3A. This result is largely driven by the combinations 
of m and g that minimize the phenotypic variance in rela- 
tively stable environments (Fig. 4). In this regard, the 
grandparental effects operate in a parameter space deter- 
mined by the same processes as the parental effects and 
therefore feed into the same variance minimizing process. 
When the parental effects are positive, the system then 
favors antagonistic, even more negative grandparental 
effects to provide the counterbalance and constrain the 
phenotypic variance. 

For simplicity, we assume that both m and g axe con- 
stant and not context dependent (Rossiter 1996; Plaistow 
and Benton 2009). For example, Plaistow et al. (2006) 
showed how the persistence (and effect size) of great- 
grandmaternal, grandmaternal, and maternal effects dif- 
fered between high- and low-food environments in the soQ 
mite Sancassania berlesei. Although work is ongoing to 
study evolvable context dependence, it will, given the com- 
plexity of the expressions derived herein assuming fixed m 
and g (equations 7 & 8), prove particularly challenging to 
incorporate multigenerational effects in this framework. 
We nevertheless built on the existing quantitative genetic 
architecture because it neatly incorporates within- and 
transgenerational plasticity. Understanding how organisms 
flexibly adjust their phenotypes to match their environ- 
ment is best achieved through concurrent investigation of 
both within-generational phenotypic plasticity and trans- 
generational plasticity via maternal and multigenerational 
effects (Uller 2008; Ezard et al. 2014). Although our inter- 
est here is how maternal and grandmaternal phenotypes 
influence the same phenotype in the focal generation (e.g.. 
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Falconer 1965, for the case of maternal effects only), the 
overarching concept that transgenerational effects alter 
phenotypic variances extends into multiple dimensions 
(Townley and Ezard 2013; Kuijper et al. 2014). The role of 
dynamic phenotypic variance is key to understanding the 
multigenerational contributions of mothers and grand- 
mothers to phenotypic evolution (Figs. 3, 4). 
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